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ABSTRACT 

A  packet  switching  network  has  the  desirable  feature  of 
rapidly  handling  short  (bursty)  messages  of  the  type  often 
found  in  computer  communication  systems.   In  evaluating  packet 
switching  networks,  the  average  time  delay  per  packet  is  one 
of  the  most  important  measures  of  performance. 

The  problem  of  message  routing  to  minimize  time  delay  is 
analyzed  here  using  two  approaches,  called  "successive  satura- 
tion" and  "max-slack" ,  for  various  traffic  requirement  matri- 
cies  and  networks  with  fixed  topology  and  link  capacities. 
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I.   INTRODUCTION 

A.   THE  PACKET  SWITCHED  NETWORK  CONCEPT 

A  new  technique  for  data  communications  that  has  evolved 
over  10  years  is  called  PACKET  SWITCHING.   In  general,  commu- 
nication networks  may  be  conveniently  divided  into  three 
types:   CIRCUIT  SWITCHING,  MESSAGE  SWITCHING  and  PACKET 
SWITCHING.   Both  message  and  packet  switching  uses  a  tech- 
nique known  as  store  and  forward  transmission. 

A  circuit  switching  network  provides  service  by  setting 
up  a  total  path  of  connected  lines  from  the  source  to  the 
destination  of  the  call.   This  complete  circuit  is  set  up  by 
a  special  signaling  message  that  threads  its  way  through  the 
network,  seizing  channels  in  the  path  as  it  proceeds.   After 
the  path  is  established,  a  return  signal  informs  the  source 
that  data  transmission  may  proceed,  and  all  channels  in  the 
path  are  then  used  simultaneously. 

The  entire  path  remains  allocated  to  the  transmission, 
and  only  when  the  source  release  the  circuit  will  all  these 
channels  be  returned  to  the  available  pool  for  use  in  other 
paths.   Circuit  switching  is  the  common  method  for  telephone 
systems. 

In  message  switching,  only  one  channel  is  used  at  a  time 
for  a  given  transmission.   The  message  first  travels  from 
its  source  node  to  the  next  node  in  its  path,  and  when  the 
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entire  message  is  received  at  this  node,  then  the  next  step 
in  its  journey  is  selected.   If  this  selected  channel  (link) 
is  busy,  the  message  waits  in  a  queue,  and  finally  when  the 
channel  becomes  free,  transmission  begins.   Thus,  the  message 
"hops"  from  node  to  node  through  the  network  using  only  one 
channel  at  a  time,  possibly  queueing  at  busy  channels,  as  it 
is  successively  stored  and  forwarded  through  the  network. 

Packet  switching  is  basically  the  same  as  message  switch- 
ing except  that  the  messages  are  decomposed  into  small  equal 
pieces  called  packets,  each  of  which  has  a  minimum  length. 
These  packets  are  numbered  and  addressed  and  make  their  way 
through  the  net  independently  of  each  other.   Thus,  many 
packets  of  the  same  message  may  be  in  tranmission  simulta- 
neously, giving  one  of  the  main  advantages  of  packet  switching 

With  packet  switching  systems,  information  is  exchanged 
in  the  form  of  short  packets.   A  packet-switched  network  can 
handle  several  different  types  of  traffic  concurrently.   These 
include  HIGH-THROUGHPUT  traffic,  for  example,  the  transmission 
of  large  data  files  between  computers,  for  which  accuracy  and 
high  average  data  speed  are  the  most  important  performance 
requirements;  LOW-DELAY  traffic,  for  example,  interactive 
communication  between  a  person  at  a  terminal  and  a  remote 
computer,  for  which  accuracy  and  low  average  message  delay  are 
important;  and  REAL-TIME  traffic,  for  example,  packetized 
speech  for  which  the  performance  of  circuit-switched  connec- 
tion must  be  approached  by  maintaining  a  relatively  constant 
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data  speed,  but  for  which  extreme  accuracy  is  not  important 
owing  to  the  redundancy  of  the  information. 

The  packet  switched  network  is  designed  primarily  for 
computer  to  computer  communication.   It  has  a  much  more 
rapid  response  which  matches  the  internal  behavior  of  compu- 
ters and  handles  information  in  much  the  same  way  as  does  a 
computer.   At  the  same  time  it  can  readily  match  the  speed 
of  attached  computers  to  that  of  the  terminal  users,  by  virtue 
of  its  internal  storage. 

The  prime  purpose  of  store  and  forward  packet  switching 
is  to  enable  communications  resources  to  be  used  effectively 
and  in  such  a  way  that  they  may  be  shared  by  many  users  opera- 
ting in  an  intermittent  fashion,  giving  each  user  a  rapid 
response  from  the  communication  network  just  at  the  instant 
when  this  is  required. 

If  there  is  need  for  transmitting  a  long  continuous  stream 
of  data,  then  a  circuit  switched  connection  makes  good  sense. 
On  the  other  hand,  if  the  data  flow  is  bursty,  then  some  form 
of  resource  sharing  can  be  used  to  great  advantage;  packet 
switching  is  an  effective  choice  here. 

Since  packets  are  stored  as  they  pass  through  switching 
nodes,  it  is  possible  to  conduct  speed,  format  and  code  con- 
version during  the  switching  process.   Another  feature  of 
packet  switching  is  its  ability  to  adaptively  select  good 
paths  for  packet  journeys  as  a  function  of  the  network 
congestion. 
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Besides  providing  small  network  delays,  packet  switching 
has  the  desirable  feature  of  rapidly  handling  small  messages 
in  spite  of  the  presence  of  long  messages  that  may  be  in 
transport  at  the  same  time,  this  is  because  of  the  decomposi- 
tion of  long  messages  into  packets.   Another  useful  property 
of  this  decomposition  is  that  the  nodal  storage  requirement 
is  reduced. 

In  evaluating  packet  switching  networks,  the  delay, 
throughput,  cost  and  reliability  are  important  measures. 

Theoretical  studies  have  been  directed  to  the  queueing  and 
network  flow  problems  in  general  and  more  specifically  to  such 
problems  as  delay  analysis,  route  assignment,  topological 
design  and  flow  control,  etc. 

The  topic  chosen  in  this  study  is  centered  around  optimal 
flow  and  minimum  capacity  assignment  and  delay  analysis  for 
packet-switched  networks.   The  networks  under  consideration 
have  a  relatively  complicated  structure  with  a  large  number 
of  source-destination  node  pairs. 

B.   THE  ROUTING  PROBLEM 

A  message  routing  procedure  is  merely  a  decision  rule 
that  determines  the  node  a  message  will  next  visit  in  its 
path  through  the  network.   The  objective  of  the  routing 
procedures  is  to  transport  packets  on  a  minimum  delay  path 
from  source  to  destination. 

In  discussing  routing  policies  for  networks,  an  important 
distinction  must  be  made  between  static  and  adaptive  policies. 
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This  distinction  depends  on  the  environment  in  which  a  policy 
is  designed  to  operate.   If  network  topology  is  not  subject 
to  changes  (due  to  failure,  modifications,  growth)  and  traffic 
inputs  are  stationary,  then  the  optimal  routing  solution  is  a 
static  solution  consisting  of  a  set  of  fixed  paths  between  all 
node  pairs.   The  traffic  between  each  source  and  destination 
pair  may  be  distributed  on  several  paths  simultaneously  in 
well  defined  proportions,  where  the  proportions  are  fixed  in 
time . 

In  a  real  network  environment,  however,  the  topology 
changes  with  time  and  user  traffic  requirements  tend  to  fluc- 
tuate more  or  less  rapidly.   To  minimize  delay  it  is  then 
necessary  to  implement  an  adaptive  routing  policy  that  can 
react  and  adjust  to  various  changes.   The  adaptivity  of  a 
policy  can  be  measured  in  terms  of  its  response  time.   A 
reasonable  procedure  is  to  use  a  periodically  refreshed 
static  routing  solution.   Indeed,  there  is  a  continuum  of 
solutions  between  the  two  extremes,  static  and  adaptive, 
characterized  by  different  response  times  and  used  for 
different  applications. 

Beside  their  use  in  operational  networks,  static  routing 
policies  find  an  important  application  in  the  network  design 
process.  During  this  process,  for  a  given  traffic  pattern  a 
prediction  of  the  throughput  and  delay  performance  of  a 
given  topolocy  is  needed.  The  routing  policy  clearly  has  a 
major  impact  on  such  performance.  Most  routing  implementa- 
tions are  adaptive,  and  unfortunately  the  analysis  of  adaptive 
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routing  policies  is  an  extremely  difficult  task.   To  simplify 
the  problem,  the  traffic  pattern  is  usually  approximated  with 
a  stationary  pattern,  and  the  routing  policy  with  a  static 
routing  policy. 

1 .   Static  Routing 

a.   Routing  Table  Representation 

A  static  routing  policy  is  represented  by  a  set 
of  routing  tables,  one  for  each  node,  indicating  how  the 
traffic  arriving  at  that  node  must  be  routed  on  the  outgoing 
links,  depending  on  its  final  destination.   The  routing  table 
for  node  i  is  an  Nx  L  matrix  Pi(n,k) ,  where  N  is  the  number 
of  nodes  in  the  network  and  Pi(n,k)  is  the  fraction  of  traffic 
directed  to  node  k  which,  upon  arrival  at  node  i,  is  routed 
through  neighbor  node  n  (Fig.  1.1).   By  the  definition, 

Pi  (j,k)  =  1 
Where  n  is  the  number  of  neighbors  of  node  i. 

The  actual  distribution  of  the  incoming  traffic 
to  outgoing  links  may  be  done  randomly  using  as  probability 
weights  the  values  Pi(j,k) .   If  for  any  node  i  there  is  only 
one  permissible  outgoing  link  from  i  to  any  destination  k,  then 

Pi(j,k)  =    1   for  neighbor  node  j. 
0   otherwise 

Such  a  routing  policy  is  referred  to  as  a  single  path 
routing  policy,  since  only  one  path  is  used  from  any  node  i 
to  any  destination  node  k.   In  general,  the  optimal  static 
routing  solution  is  a  multipath  solution,  allowing  for  the 
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simultaneous  use  of  several  routes  in  order  to  minimize 
delays.   The  routing  table  representation  may  also  be  used 
for  adaptive  policies. 


•o 


l 
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Pi(l,k) 
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. ■■ j 
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Figure  1.1.   Routing  Table. 

b.   Problem  Formulation 

In  packet  switched  network,  messages  are  seg- 
mented into  small  packets,  and  each  packet  travels  from  its 
source  to  its  destination  via  a  set  of  intermediate  nodes. 
While  awaiting  transmission  to  the  next  node  in  the  route, 
the  packet  must  be  stored  until  the  link  to  that  node  is 
free.   Thus,  at  each  node  there  are  several  queues,  one  for 
each  output  channel. 
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Packet  flow  requirements  between  nodes  arise  at 
random  times;  therefore,  link  flows,  queue  length,  and  packet 
delay  are  random  variables.   A  static  routing  problem  can  be 
defined  as  the  problem  of  finding  the  static  routing  policy 
which  optimizes  the  flow  and  minimizes  the  average  time  delay. 

Considering  the  relationship  between  routing  policy 
and  link  flows,  the  problem  can  be  formulated  for  solution  of 
optimal  flow  and  minimum  capacity  with  respect  to  the  given 
requirement  as  follows, 
given:   (a)  topology 

(b)  requirement  matrix  R 
minimize:       alpha 
subject  to:   (a)  conservation  of  flow 

(b)  flow   capacity  x  alpha 
Where  alpha  is  a  nonnegative  scale  factor. 

The  mathematical  model  and  solution  techniques  are 
discussed  in  detail  later. 
2 .   Adaptive  Routing 

The  adaptive  routing  problem  consists  of  defining  a 
procedure  that  dynamically  updates  routing  tables  according 
to  changes  in  the  network.  This  procedure  will  include  the 
acquisition  of  the  status  of  all  neighboring  nodes,  which  a 
packet  associated  with  the  source  and  destination  node  pair 
could  visit  next  after  being  processed  at  the  current  node. 
The  neighboring  nodes  are  those  with  direct  channels  linking 
to  a  given  node. 
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The  structure  of  the  routing  table  and  the  ways  the 
information  contained  therein  are  to  be  used  vary  according 
to  the  individual  routing  procedure.   For  example,  if  single 
path  routing  is  used,  only  one  path  links  each  source  and 
destination  node  pair;  therefore,  only  a  single  neighboring 
node  associated  with  each  node  pair  is  listed  in  the  routing 
table . 

The  design  of  routing  is  a  process  of  determining  the 
structure  of  routing  tables  and  specifying  the  procedures  of 
using  them.   The  choice  of  routing  depends  on  many  factors 
such  as  network  topology,  available  facilities,  throughput 
and  delay  requirements,  and  so  on.   Adaptive  routing  tech- 
niques are  ruled  out  in  this  study. 
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II.   OPTIMAL  FLOW  AND  MINIMUM  CAPACITY  ASSIGNMENT 

The  optimal  flow  is  the  best  selection  of  paths  between 
source  and  destination,  given  the  traffic  requirements  and 
the  network  configuration.   The  definition  of  best  paths  may 
vary  depending  on  the  nature  of  the  traffic.   In  this  study, 
the  best  paths  will  be  regarded  as  the  paths  of  minimum 
"distance",  where  distance  can  be  interpreted  in  various  ways, 
for  instance  as  time  delay  along  the  links. 

To  investigate  average  time  delay  of  the  networks,  two 
different  approaches  are  discussed  using  classical  optimiza- 
tion procedures . 

A.   NETWORK  MODELING 

In  a  distributed  packet-switched  network,  packets  are 
transmitted  over  a  collection  of  switching  nodes  in  tandem. 
In  modeling  a  packet  switched  network,  the  objective  is 
usually  to  optimize  the  total  flow  of  the  networks,  to 
minimize  the  capacity  of  each  link,  or  to  maximize  the 
utilization  of  the  network.   The  problem  may  be  formulated 
as  follows. 

given  topology 

requirement  matrix  R 

capacity  of  each  link  (arbitrary) 

minimize:         alpha 

subject  to:       (1)  conservation  constraints. 
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lFl(k)  -    lFl(h)  =        Ri(k) 

outgoing       incoming        input 

links  links        at  node  i 

from  node  i    to  node  i       for  node  k 

(2)  capacity  constraints 
Ft (k)    alpha  x  CC 

(3)  positive  constraints 
iFl (k)  >  0 

Where  F£(k)  is  the  flow  on  link  destined  for  node  k.   We 

next  look  at  each  constraint  further  in  detail. 

1 .   Conservation  Constraints 

At  each  node  i,  the  total  outgoing  packet  flows 

destined  for  node  k  equal  the  incoming  flows  plus  the  input 

to  be  transmitted  to  node  k. 

iFlCk)  -        lFl(k)  =        Ri(k) 

outgoing       incoming        input 

links  links        at  node  i 

from  node  i    to  node  i       for  node  k 

where 

k  =  destination  node  (k  =  1,2,3...N). 

N  =  number  of  nodes. 

L  =  number  of  links. 

i  =  source  node  (i  =  k,  i  =  1,2,3. . .N). 

i   -    link  between  two  nodes  (t =  i,2,3...L). 

FC(k)  =  flow  on  link  destined  for  node  k. 

Ri(k)  =  input  destined  for  node  k  at  node  i. 
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2 .  Capacity  Constraints 

The  total  flow  destined  for  node  k  on  a  certain  link 
should  be  less  than  or  equal  to  the  capacity  of  the  link  t. 
Since  linear  programming  algorithms  have  difficulties  dealing 
with  inequalities,  slack  variables  are  added  to  the  constrains 
to  absorb  unused  resources  and  thus  to  force  equalities  to 
appear. 

When  slack  variables  are  introduced,  the  problem  in 
standard  form  becomes. 

^  FC(k)  -  alpha  X  Cl  +   SI   =  0 
k=l 
where 

Cl   =   capacity  of  the  link. 

alpha  =  constant  to  optimally  minimize  the  link 

capacity  and  total  flow  on  link. 
SZ   =  slack  variables. 

3 .  Positive  Constraints 

If  we  consider  the  network  to  have  unidirectional 
links,  the  bidirectional  flows  between  two  nodes  are  divided 
into  two  one-way  links  in  opposite  direction,  and  the  flow  on 
each  link  is  positive. 

Ft  (k)  >  0 

Using  matrix  notation,  the  model  simply  may  be 
expressed  as  a  linear  programming  problem. 

minimize:      alpha 

subject  to:    AX  =  R 
x   ^>  0 
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B.   MODEL  PROGRAMMING 

In  order  to  use  a  readily  available  Mathematical  Program- 
ming System  such  as  MPSIII,  which  can  solve  linear  programming 
problems  with  up  to  4000  rows  and  theoretically  with  an  un- 
limited number  of  variables,  the  major  programming  problem 
with  large  scale  networks  is  how  to  generate  the  input  data 
for  the  MPSIII  program.   In  other  words,  how  do  we  generate 
the  matrix  form  of  the  model. 

1.   Node  to  Link  Incidence  Matrix 

For  any  network  processing  to  take  place  in  the  compu- 
ter, the  network  structure  must  be  represented  in  some  machine 
understandable  form.   A  wide  spectrum  of  different  represen- 
tation methods  are  available.   Most  of  these  are  based  on  an 
incidence  matrix  scheme.   A  matrix  form  is  convenient  for  the 
analysis  of  networks  since  the  resulting  matrices  are  in  a 
format  suitable  for  mathematical  analysis. 

The  network  of  Fig.  2.1  will  be  used  to  illustrate  the 
method.   The  node  to  link  incidence  matrix  E(i,£)  describes 
the  links  connected  to  the  node  such  that 

E(i,£)  =f  1   if  link  I    is  outgoing  from  node  i. 
-1   if  link  t    is  incoming  to  node  i. 
„  0   if  link  t    is  not  connected  to  node  i. 
The  incidence  matrix  for  input  at  node  i  destined  for  node  k 
is  a  modification  of  the  node  to  link  incidence  matrix  E(i,£) 
according  to  the  following  rules: 
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Ek(i,£)  =  [  0    for  columns  corresponding  to  links 

outgoing  from  node  k. 
delete  the  row  if  i  =  k 
E(i,.£)  otherwise 
For  the  illustration  network  (Fig.  2.1). 


1 

-1 

1 

-1 

0 

-1 

1 

0 

0 

-1 

0 

0 

-1 

1 

1 

for  k  =  1 


1 
0 


-1 

1 


E2(i,£)  and  E3(i,£)  may  be  represented  in  the  same  manner 
for  k  =  2 


-1 
0 


1    -1 

-1     1 


for  k  =  3 
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E(i,.£)  is  a  3x5  matrix  and  Ek(i,£)  is  a  2x5  matrix  in  size. 


Figure  2.1.   Illustration  Network. 


2 .   Matrix  Representation 

Matrix  representation  is  a  convenient  tool  for 
programming.   The  complete  matrix  representation  for  the 
illustration  network  is  shown  in  Appendix  A. 

The  general  model  matrix  AX  =  R  may  be  represented 
as  follows: 


El 


E2 


A 


E3 


II   12   13 


Ek 
Ik 


-  Ci 


Fl(l) 

Ft  (2) 

f£(3) 

Ft  (k) 

Ft  (k) 

SI 

Alpha 

Rid) 
Ri(2) 
Ri(3) 


Ri(k) 


01 


24 


where 

Ek     =  (N-l)  x  L  incident  matrix. 

Ik     =  L  x  L  identity  matrix  with  zero  elements 

corresponding  to  zero  columns  of  the  Ek  matrix. 
I      =  L  x  L  identity  matrix. 
Ri (k)  =  requirement  matrix  which  is  the  inputs  to  be 

sent  from  source  node  i  to  destination  node  k. 
The  size  of  the  model  in  matrix  representation  is; 
A  =  ((N-l)  x  N  +  L)  x  ((N+l)  x  L  +  1) 
X  =  (  (N+l)  x  L  +  1)  x  1 
R  =  (  (N-l)  x  N  +  L)  x  1 
For  example,  the  network  consisting  of  20  nodes  and 
40  links  has  the  matrix  which  is  A  =  420  x  841,  X  =  341  x  1 
and  R  =  420  x  1  in  size. 
3.   Programs 

To  generate  the  MPSIII  input  data  and  to  solve  the 
model,  requires  two  programs  which  are  listed  at  the  end  of 
this  study  (Appendix  C  and  D) . 

The  first  one  is  the  data  (model)  generation  program 
(Appendix  C)  written  in  the  Fortran  language.   The  output  of 
this  program  uses  special  notation  to  designate  the  flow  of 
the  link  connecting  node  pairs,  destination  node  of  the  flow, 
flow  variables,  slack  variables  and  input  at  the  node.   These 
notations  are  composed  of  a  letter  which  is  one  of  C,  L,  X  or 
S  and  7  numerical  digits. 
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The  notation  Ci,  LZ,    XZ,    and  St  imply  conservation 
constraints,  capacity  constraints,  flow  variables,  and  slack 
variables  respectively.   The  rest  of  the  six  numbers  following 
CZ,    LZ ,    XZ   and  SZ    indicates  link  connecting  node  pair  with 
first  4  digits  and  destination  node  or  slack  variable  number 
with  last  2  digits.   For  example,  the  notation  C1000201  indi- 
cates the  conservation  constraint  row  which  is  the  input  at 
node  2  destined  for  node  1,  L1020100  the  link  capacity  row  for 
the  link  connecting  between  node  2  and  node  1,  X1020101  the 
flow  variable  on  the  link  connecting  between  node  2  and  node 
1  destined  for  node  1  and  S1000001  the  first  slack  variable. 

The  second  program  is  the  MPSIII  program  developed  by 
Management  Science  Systems,  Inc.   It  provides  elaborate  con- 
trol language  for  the  formulation  of  solution  strategies  for 
mathematical  programming  problems. 

Individual  instructions  in  this  language  bring  in 
quite  elaborate  sections  of  code  designed  to  execute  the  step 
as  efficiently  as  possible. 

Using  the  special  notation  discussed  above,  there  is 
a  limitation  on  the  size  of  the  problem  which  can  be  solved, 
because  data  names  must  consist  of  8  characters  maximum. 
This  means  that  the  maximum  network  size  must  be  less  than 
100  nodes,  if  the  computer  capacity  allows. 

C.   SOLUTION  TECHNIQUES 

As  mentioned  before  ,  the  classicial  optimization  proce- 
dures to  solve  the  investigation  model  provide  solutions 
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which  are  the  optimal  paths  and  minimum  link  capacities  for 
the  network.   Different  optimal  solutions  may  be  expected 
depending  on  different  methods  of  optimization.   The  two 
different  approaches,  each  involving  iteration  of  the  linear 
programming  solution  procedure  are  discussed  below,  and  illus- 
trated in  terms  of  one  example  network  problem.   These 
approaches  are  based  on  the  characteristics  of  linear  pro- 
gramming problems. 

1.   Successive  Saturation  Approach 

The  standard  procedure  of  the  MPSIII  program  listed 
at  the  end  of  this  study  (Appendix  D  and  E)  produces  an 
"optimal"  solution  for  alpha  from  the  output  data  of  the  data 
generation  program  in  the  first  run.   But  the  corresponding 
flow  solution  is  not  a  good  solution  in  terms  of  the  average 
time  delay  in  the  packet-switched  network,  when  it  is  ana- 
lyzed by  the  Kleinrock's  delay  analysis  model  (Ref .  4) .   The 
problem  is  that  only  the  flow  thru  the  saturated  links  is 
optimum. 

In  order  to  reduce  the  time  delay,  Iteration  runs  are 
attempted  until  all  links  are  saturated.   The  detailed  itera- 
tion method  for  the  example  network  (Fig.  3.1)  which  has  5 
nodes  and  6  links,  each  with  original  link  capacities  10  (to 
simplify  the  problem)  is  discussed  below. 


27 


Figure  3.1.   Example  Network. 

The  ouput  of  MPSIII  program  is  divided  into  row 
section  and  column  section  (Ref .  9  and  Appendix  D  and  E) . 
The  objective  function  value  (which  is  the  same  as  the  alpha 
value)  and  the  links  saturated  with  the  alpha  value  in  the 
capacity  constraint  rows  indicated  with  the  first  letter 
notation  L  can  be  read,  and  the  flows  on  each  link  and  the 
alpha  value  in  column  section.   In  the  example  network,  for 
the  first  iteration  the  saturation  level  (alpha  value)  is 
0.25,  and  the  saturated  links  are  L1010300  (link  3)  and 
L1040300  (link  7).   Their  activity  and  slack  activity  levels 
are  all  zero,  but  their  dual  activity  level  is  not  (pp  48) . 

For  the  next  iteration,  the  changed  input  datas  are 
that  the  original  capacity  multiplied  by  the  saturation  level 
of  first  run  are  moved  from  the  alpha  column  to  the  right 
hand  side  (RHS)  column.   In  the  example,  the  new  link  capaci- 
ties 2.5  (0.25  x  10)  are  moved  to  RHS  (pp.  50  ),  and  the 
problem  is  solved  again.   Repeating  this  iteration  procedures 
until  all  links  are  saturated  and  moved  to  RHS,  the  objective 
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function  value  becomes  zero.   When  the  objective  function 
value  is  zero,  the  final  flow  solution  is  the  desired  solution 
of  this  successive  saturation  approach.   In  the  example  case, 
the  six  iterations  are  required  to  solve  the  problem.   These 
whole  procedures  except  programs  and  input  datas  from  2nd  run 
for  the  example  network  are  listed  at  the  end  of  this  study 
(Appendix  D) . 

2 .   Max- Slack  Approach 

To  simplify  the  solution  of  the  above  successive  satu- 
ration approach,  another  approach  called  the  max-slack  approach 
is  considered.   It  requires  only  two  iteration  of  the  program, 
but  is  otherwise  quite  similar.   Like  in  the  successive  satu- 
ration approach,  after  the  first  run,  all  saturated  link 
capacities  are  multiplied  by  the  alpha  value  of  the  first  run 
are  moved  to  the  RHS,  and  in  addition,  all  slack  variable 
columns  are  summed  as  the  objective  function  in  the  data  card 
for  the  second  run. 

In  the  example  network,  eight  slack  column  variables 
are  added  on  the  objective  rows,  and  all  new  capacities  2.5 
(0.25  x  10)  are  moved  to  RHS.   This  is  clearly  shown  on  the 
output  of  control  language  "picture",  and  the  program,  input 
data  and  the  output  are  also  listed  at  the  end  (Appendix  E) . 

D.   DELAY  ANALYSIS 

The  average  time  delay  of  a  network  is  the  average  time 
a  packet  spends  in  the  network  traveling  from  its  source  to 
its  destination  nodes.   To  obtain  a  tractable  expression  for 
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average  time  delay  T  in  terms  of  the  capacity  of  link  I,    it 
is  necessary  to  make  simplifying  assumptions.   Packet  lengths 
are  assumed  to  be  exponentially  distributed,  and  are  re-chosen 
with  statistical  independence  at  each  node  to  form  a  Poisson 
process . 

These  assumptions  known  as  the  'Independence  assumption', 
introduced  by  D.  Kleinrock  (Ref.  4),  make  each  queuing  prob- 
lem independent.   In  the  analysis,  each  queue  in  a  packet 
network  is  assumed  to  be  an  M/M/l  system  with  arrival  rate 
of  Gt    and  to  be  independent  of  each  other. 

The  average  delay  for  link  I    is  given  by  the  waiting  time  in 
the  queue  as 

T   =     1 

ll  uC£-Ge 

Where  Gt    is  the  packet  traffic  in  the  link  and  1/u  the 

average  packet  length  in  bits.   To  form  a  suitable  average 

over  all  the  queuing  processes,  the  T   are  weighed  by  Gi/r 

where  r  is  the  total  input  packet  rate  to  the  network.   In 

this  way,  the  total  delay  suffered  by  all  packets  per  second 

of  network  operation  I    Gl-T    is  divided  bv  the  total  number 

I 
of  packets  carried  by  the  network  per  second. 

The  total  average  time  delay  per  packet  through  the  net- 
work is  then  given  by 

_  i    i       Gi       __  i  i       Fi 


r  t=i  uCr  Gi     rt=i  ct  "  Ft 
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Where  Gt/u  =   Ft,    and  L  is  the  total  number  of  links  in 
the  system  and  1/u  is  equal  to  the  average  flow  on  link  c. 
It  should  be  noted  that  the  value  of  C£  (the  minimum  required 
capacity)  must  be  greater  than  the  total  of  all  the  flows  on 
link  t.       In  other  words,  the  value  of  C£  must  be  greater  than 
the  value  which  is  the  first  run  saturation  level  multiplied 
by  the  arbitrary  capacity  C  for  the  first  run. 

The  value  of  T  includes  both  of  the  waiting  time  and  the 
service  time.   The  waiting  time  is  subject  to  the  interference 
of  all  other  traffic  within  the  network  that  consists  of  the 
data  traffic  as  well  as  the  control  traffic.   This  average 
time  delay  can  be  used  to  analyze  the  relative  goodness  of  the 
routing  strategies  considered  in  this  study. 

In  deriving  this  basic  formula,  a  number  of  factors  (such 
as.  processing  time  and  propagation  delay)  are  neglected.   In 
any  realistic  network,  these  variables  as  well  as  others  must 
be  considered  in  the  analysis  of  the  networks.   According  to 
the  above  formula,  the  average  time  delay  depends  only  on  the 
aggregate  flow  of  the  links.   The  reduction  of  the  flow  Ft    of 
each  link  results  in  the  reduction  of  the  average  time  delay 
for  the  network  which  has  constant  input  packet  rate  and  link 
capacities. 

For  the  example  network  (Fig.  3.1),  the  arbitrary  capa- 
city is  10  and  the  value  of  Ct    must  be  greater  than  2.5 
(0.25  x  10) .   The  average  time  delay  for  the  capacity  Ft  =  10  , 
and  the  input  (requirement  matrix)  of  2  units  from  node  1  to 
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node  4  and  5  units  from  node  2  to  node  3  are  Tss  =  2.12/r 
for  the  successive  saturation  case  and  Tms  =  2.0  83/r  for  the 
max-slack  case.   To  make  these  calculations  easy,  the  data 
summarizing  the  final  output  run  are  listed  in  Table  I.   We 
see  that  there  is  not  much  advantage  of  one  procedure  over 
the  other  in  terms  of  performance,  whereas  the  max-slack 
solution  requires  only  two  iterations,  rather  than  6  for  the 
successive  saturation  algorithm. 
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TABLE  I 
LINK  FLOW  FOR  THE  EXAMPLE  NETWORK 


(    ) ;  Max-Slack 


LINK        DESTINATION  NODE TQTAL 

(Nodes)        12       3       4       5        FLOW 


1-2 

(2) 
2 

(2) 

z. 

2-1 

(2.5) 
2.5 

(2.5) 
2.5 

1-3 

(2.5) 
2.5 

(2.5) 
2.5 

3-4 

2-4 

(2.5) 
2.25 

(2.5) 
2.25 

2-5 

0.25 

(2) 
2 

(2) 
2.25 

4-3 

(2.5) 
2.5 

(2.5) 
2.5 

5-4 

0.25 

(2) 
2 

(2) 
2.25 

Req 

Mat 

C=10, 

RX(4): 

=  2,    R2 

(3)  =5 
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III.   EXPERIMENTAL  RESULTS  AND  CONCLUSIONS 

A.   EXPERIMENTS  AND  COMPARISON  OF  THE  RESULTS 

The  experiment  was  conducted  for  two  experimental  networks 
to  analyze  both  approaches,  in  addition  to  the  example  network 
These  networks  are  specified  in  Appendix  B.   The  networks  were 
used  in  the  experiment  with  requirement  matrices  of  15  units 
from  node  1  destined  for  node  5,  9  units  from  node  7  destined 
for  node  2  for  the  9  nodes/36  links  network  (Fig.  B.l) :   and 
15  units  from  node  2  destined  for  node  12,  4  units  from  node 
5  destined  for  node  6,  and  10  units  from  node  10  destined  for 
node  1  for  the  13  nodes/60  links  network  (Fig.  B.2;.   Both  the 
successive  saturation  and  max-slack  approaches  were  used  for 
the  successive  saturation  approach,  the  21  iterations  are 
required  in  this  experiment. 

The  results  and  their  comparison  is  listed  in  Tables  II, 
III,  and  IV.   These  are  obtained  by  writing  each  flow  variable 
value  from  the  final  output  run  on  each  link  of  the  networks, 
and  summing  these  flow  variable  values  to  obtain  the  total 
flow  on  the  link  (Appendix  F  and  G) . 

As  shown  in  Tables  II  and  III,  the  number  of  iterations 
required  and  the  total  link  usage  (utilization)  are  much 
greater  for  the  successive  saturation  approach,  but  fewer 
large  link  capacities  are  required  than  in  the  max-slack 
approach.   The  average  time  delay  depends  on  the  capacity 
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value  C£  (Table  IV) .   When  the  capacity  value  Ct   is  close 
to  the  minimum  required  capacity,  the  time  delay  is  less  for 
the  successive  saturation  approach  than  for  the  max-slack 
approach  in  both  networks,  and  vice  versa.   It  is  significant 
for  the  large  scale  network,  and  for  large  enough  link  capa- 
cities Ql ,    that  the  time  delay  difference  becomes  negligible. 

Also,  in  the  process  of  this  study,  it  was  noted  that  both 
approaches  yield  the  same  aggregate  flows  in  the  case  of  very 
simple  networks  (nodes  less  than  or  equal  to  4) ,  and  for  com- 
pletely symmetric  networks  and  requirement  matrices  (input  and 
capacity  matrices) .   For  example,  the  time  delays  are  the  same 
for  both  approaches  with  2  units  input  from  the  outmost  nodes 
to  the  directly  opposite  nodes  respectively  (1  to  13,  2  to  12, 
5  to  9 ,  12  to  2  and  13  to  1)  for  the  13  node  network,  because 
all  links  total  flow  are  3/4 .   The  detail  experimental  results 
are  listed  in  Appendix  H. 

B.   CONCLUSIONS 

This  is  an  elementary  study  of  the  routing  design  problem 
for  a  packet  network:   given  a  traffic  requirement  matrix, 
minimize  the  average  time  delay  per  packet,  subject  to  finding 
a  feasible  flow  for  a  network  with  fixed  topology  and  link 
capacities . 

The  general  performance  characteristics  and  advantages 
of  the  two  approaches  are  investigated  for  small  and  simple 
networks.   This  illustrates  one  way  to  solve  the  design 
oroblem  under  limited  conditions. 


37 


TABLE  IV 
TIME  DELAY  FOR  THE  CAPACITIES 


CAPACITY       4        4.2      4.5      5       10       15 

SUCC  244. 1/r  109. 5/r  59.9/r  10 . 74/r  5 . 82/r 

9/36   SAT 

Net    ^X 

28.56/r  116. 6/r  59. 5/r  10.41/r  5.73/r 

SLACK 

SUCC    169. 9/r   115 . 6/r  15. 1/r      60. 1/r  15 . 34/r  8.93/r 

13/60  SAT 

MAX 
Wetl     -      335.  7/r   188.  3/r  114. 2/r  69.3/r  14.29/r  7.99/r 
SLACK 
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APPENDIX  A 
MATRIX  REPRESENTATION  FOR  ILLUSTRATION  NETWORK 
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Figure  A.l.   Matrix  Representation  for  Illustration  Network 
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APPENDIX  B 
EXPERIMENTAL  NETWORKS 


Figure  B.l.   The  9  Node/36  Link  Network 


Figure  B 


.2.   The  13  Node/60  Link  Network, 
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APPENDIX  C 
MODEL  GENERATION  PROGRAM  FOR  EXAMPLE  PROBLEM 

$J03  1538P 

CSOPTIQNS    FREE 

:    THIS    IS    THE     DATA    GENERATION    PROGRAM    FOR    MPS! i I. 

C    THE    FORMAT    SHGJLD    BE    CHANGED    ACCORDINGLY. 

C    N    IS    NUMBER    OF    NOCES,     L    IS    NUMBER    OF    LINKS. 

C    ItdiU     IS    OEFINEO  AS    FOLLOWS: 

C  l.LINK    INCOMING    TO    NOD  E I    IS   -1. 

C         2.  LINK    OUTGGING   FROM   NODE1    IS     1. 

C         3.  LINK    UNCONNECTED     WITH   NODEI     IS    0. 

C    IZCOLU,L)     IS    DEFINED    AS     FOLLOwS: 

C         l.LINK    OUTGOING    FROM   NODE    K    IS    0. 

C         2. OTHERWISE,    1. 

C    IRUtKJ     IS    THE    INPUT    MATRIX    FROM    NODE    I    TO   NODE    K 

C    CAP(L,1)     IS     ThE    LINK     CAPACITY    MATRIX. 

C    COLNKltJl)     IS    THE  LINK    MATRICX    CONNECTING    NODES. 

r 
U 

C 

DIMENSION    INC  (13,  60)  , IE (13,6 0), I  ZC0L(13,  oO)  , 

*  IRN(13,60) , IR( 13, 13  J, CAP (6 0,1  J  ,  ICOLN1 1  1 ,60 ) 

*  ICOL,N(  13,60), IC0L\I2(13, 60) 

INTEGER     I  ,J,K,L,N,NJLL,  IDEM  ,  KA,  J  1 ,  INC  2, 1  A  , 

*  ID,IQ,IROWNG,IL,KC  ,KD,IRN1 
C 

C    ***    DATA    INITIALIZATION    ***** 
READ  (5,  20C)    N  ,L 

READ(5,225)(IC0LN1(1  ,J1),  J1  =  1,L) 
READ (5, 2  2  0) (( lei  I, J) ,1=1, N) , J=I,L) 
READ (5, 230) ((  IR( I  ,K)  ,  1= 1  ,N)  , K=l , N) 
READ (5, 2  40) (CAP(  IQ,1  ),IQ=1,L) 
200    F0RMAK215I 
220     F0RMAT(5I3) 
225    FORMAT(  4-110) 
230    FORMAT!  513) 
2<tO    FQRMAT(4E12.5  ) 
C 

C    ***    ZERO    COLUMN    GENERATION    **** 
DO    19    J    ■    1,L 
DO    2  0    I     =    1,N 
INCE    =    IE(I,J ) 
IFdNCE     .EQ.     1)    GO   TO    21 
INCE    =    1 
GO    TO   23 
21     INCE    =    0 
23    IZCOLU  ,  J)  =  INCE 
2  0    CONTINUE 
19    CONTINUE 
C    ***    COLJMN    NUMBER    GENERATION    **** 
IAD    =    1000000 
DO    1    K    =     1,N 
IROWNO    =    0 
DO    12    Jl    =    1,  L 

IC0LN(K,J1)    =     IC0LN1(1,J1)    +    K    *     IAD 
1CGLN2(K,J1)     =    IC0LN(K,J1)    -    K 
12    CONTINUE 
C    ***     INCIDENT    MATRICX     GENERATION    *** 
DO    2     I    =    1,N 
IF(I     .EQ.    K)     GO    TO    2 
DO    3     J    =    1,L 
NULL    =    IZCOL(  K,J) 
IF(I     .EQ.    K)     NULL    =    0 
IDEN    =    IE(I,J ) 
INC(  I*J  )     =    NULL*IDEN 
3     CONTINUE 
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2 


61 
1 


7 

C  *** 


63 


62 

65 


66 


67 

C  *** 


63 
69 


140 
145 
500 
510 
520 
540 
550 
560 
565 
570 
580 


IRGWNO 
IRN< ItK 
CONTINU 
ROW  SEC 
ID  =  1 
DO  61  I 
IF(I  .E 
wRITE(6 
CONTINU 
CONTINU 
CAPACIT 
DO  7  J 
WRITE(6 
CONTINU 
COLUMN 
WRITE*6 
DO  65  K 
DO  62  J 
NA  =  0 
DO  63  I 
IFil  .£ 
NULL  = 
I  DEN  = 
IHCi I  ,J 
INC2  ■ 
IF(INC2 
NA  ■  NA 
WRITE(6 
CONTINU 
IFCNA  . 
WRITEio 
CONTINU 
CONTINU 
CAPACIT 
DO  66  J 
IA  =  J 
WRITE(6 
CONTINU 
ALPHA  C 
WRITE  (6 
DO  67  J 
WRITE(6 
CONTINU 
RHS  GEN 
WRITE<6 
DO  6  9  K 
DO  68  I 
IF(I  .E 
WRITE16 
CONTINU 
CONTINU 
DO  5  IQ 
CONTINU 
WRITE(6 
FORM  AT ( 
FORMAT! 
FORM  AT ( 
FORM  AT ( 
FORMAT( 
FORM  AT ( 
FORMAT* 
FORM  AT ( 
FORM  AT  ( 
FORM AT ( 
FORM  AT  ( 
STOP 
END 


-     1*100 

)    =    IRCWNO    +    K    +    IAD 

TICN    GENERATION    *** 


=    1,N 
Q.    ,<) 
#140)1 

E 
E 

Y  ROW 
*  lfL 
,145) 
E 

SECTIO 
,500) 

=    1,N 
=    1,L 

=    1,N 
Q.    K) 
IZCOL( 
IECIt J 

)     =    NU 
INCd, 

•  EQ. 

♦  1 
,510)1 

EQ.    0) 
,520)1 

E 
E 

Y  GENE 
=    1,L 

+    1000 

,540) 

E 

OLUMN 

,565) 

=  1,L 
t  570) 
E 

ERATIO 
,550) 

=    1,N 

=    1,N 
Q.    K) 
,  560) 
E 
E 

=  1,L 
E 

,580) 
IX, «E 
IX, «E 
' COLUM 
4X,IX* 
4X,,XI 
4X,«SI 
•RHS* ) 
4X,'  INPUT 
4X,' ALPHA 
4 X,1 ALPHA 
■  ENOATA*  ) 


GO    TO   6  1 
RN( I  ,K) 


GENERATION    *** 

ICGLN2(  1,J) 

N  GENERATION    *** 


GO    TO   6  3 

K,J) 

) 

LL#ID£S 

J) 

0)    GO   TO    63 

COLNUi  J  ),  I  RISK  I  ,K)  ,  INC2 

GO    TO    62 
COLN(K,  J  ),1CQLi\J2U,  J),  ID 


RATION    *** 

000 
IA,ICCLN2(1,J)  ,ID 

GENERATION    *** 
ID 

1C0LN2(  1  ,J)  ,CAP(J,1) 

N   *** 


GO    TO   6  8 
IRNUt*  IflRUff&l 


C  ,  17) 
L"  ,17) 

NS' 

,17 

,17 

,17 


2X, 'C1  ,I7,2X,I5) 
2X,  •L1  ,I7,2X,I5) 
2X,  'L1  ,I7,2X,I5) 

,5X , 'C* ,17, 2X,I 5) 

,5X  ,  •OBJ1 , 7X,15 ) 

,5X ,'L* ,I7,2X,512.5) 
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SENTRY 

5 

3 

010200 

0201 

,00 

010300 

030200 

02  0400 

020500 

0^030  0 

050^00 

1    -1 

0      0 

0 

-i       1 

0       0 

0 

1       0 

-1       0 

0 

0    -1 

1      0 

0 

0       1 

0   -1 

0 

0      1 

0      0 

-1 

0       0 

-1     1 

0 

0      0 

0   -1 

1 

0      0 

0      0 

0 

0      0 

0      0 

0 

0      4 

0      0 

0 

4      0 

0      0 

0 

0       0 

0      0 

0 

-10. 

■10. 

1 

-10. 

-10. 

-10. 

•10. 

1 

-10. 

-10. 

SENTRY 
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APPENDIX  D 

OPTIMIZATION  PROGRAM,  INPUT  DATA 
MATRIX  PICTURE  AND  OUTPUT  FOR  EXAMPLE  PROBLEM 

//CHANG  JOB   (1538? 18081 #» CHANG  JOON  WGOMG ' ,Cl ASS  =  4 

//*       APSlil    EXAMPLE  PROBLEM 

//* 

//       EXEC       MSSMPS 
//CPC.SYSIN      00      * 

PROGRAM    (•  N0«  I 
INITIALZ 

TITLE     ('EXAMPLE    PROBLEM,     1ST    RUN    FOR    SUCC-SATM 
MOVE       (XQBJ.'GBJ1 ) 
MOVE        (XRHS,1  INPUT*  ) 
MOVE        (XOATA,  ■  LINEQS1  ) 
MOVE       (XPBNAMEt'PROJECT*  I 
CONVERT    <•  SUMMARY  I 
SETUP     CMINM 
BCOOUT 
PICTURE 
WHIZARO 
PRIMAL 
SOLUTICN 
EXIT 
PENO 
/* 

//EXEC.SYSIN         DO       * 
NAME  LINEQS 

ROrtS 
N  OBJ 
E  C1000201 
E  C1000301 
E  C1000401 
E  C1000501 
E  C1000102 
E  C1000302 
E  C1QQ04Q2 
E  CI  0005 02 
E  C1000103 
E  C1000203 
E  C10Q0403 
E  C1000503 
E  C1Q00104 
E  CI  0002 04 
E  C1000304 
E  C100Q504 
E  C1000105 
E  C1000205 
E  CI  0003 05 
E  C1000405 
E  L1010200 
E  L1020100 
E  L1010300 
E  L1030200 
E  L1020400 
E  L1020500 
E  LI  0403  00 
E  L1050400 
COLUMNS 

X1020101  C1J00201 
X1020x01  L1020100 
X1030201  C1000201 
X1030201  C10OO3O1 
X1030201  L1030200 
X1020401  C1J00201 
X1020401  C1000401 
X1020401       L1J2040Q 
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X1020501  C1000201 

X1020501  C1000501 

X1020501  L1020500 

X1Q40301  CI000301 

XI 0403 01  CI  000 401 

XI 0403 01  L 1040 200 

X1050401  C1300401 

X1050401  C10QQ501 

X105Q4Q1  L1050400 

X1010202  C1000102 

X1010202  L101O2OO 

X1010302  C1000102 

X1010302  C10Q02Q2 

X1010302  L1010300 

X1030202  C1000302 

X1030202  L1030200 

X1040302  C1J00302 

X1040302  C1000402 

X1040302  L1040300 

X1050402  C10Q0402 

X1050402  C1J00502 

X1050402  L1050400 

X1010203  C10001Q3 

X1010203  CI  300  203 

X1010203  L1010200 

X1020103  C13001Q3 

X1020103  C1000203 

XI 0201 03  L13201Q3 

X1010303  C1000103 

X1010303  L1010300 

X1020403  C1QQ0203 

X1020403  C1000403 

X1020403  L102Q40Q 

X1020503  C1000203 

X1020503  C1000503 

X1020503  11020500 

X1040303  C1000403 

X1040303  L1040200 

X1050403  C1000403 

XI 050403  CI  000 503 

X1050403  L1050400 

X1010204  C1000104 

X1010204  C1000204 

X1010204  L1010200 

XI 0201 04  CI  000 104 

X1020104  C1000204 

XI  0201  04  L1020100 

X1010304  C10QQ104 

X1010304  C1000304 

X1010304  L1010300 

X1030204  C1000204 

X1030204  C1000304 

X1030204  L1030200 

X1020404  C1000204 

X1020404  L1020400 

X1020504  C1000204 

X1020504  C1000504 

XI 02 05 04  L1020  500 

X1050404  C1000504 

X1050404  L1050400 

X1010205  C1000105 

X1010205  C1000205 

X1010205  L1010200 

X1020105  C1000105 

X1020105  C1000205 

X1020105  L1020100 

X1010305  C1000105 
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RHS 


XI 
XI 
XI 
XI 
XI 
XI 
XI 
XI 
XI 
XI 
XI 
XI 
XI 
SI 
SI 
Si 
SI 
SI 
SI 
SI 
SI 
AL 
AL 
AL 
AL 
AL 
AL 
AL 
AL 
AL 


IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
END ATA 


010305 

010305 

030205 

0302  05 

030205 

020405 

02040  5 

020405 

020505 

020505 

0403  05 

040305 

0403  05 

000001 

000002 

0000  03 

0000  04 

000005 

0000  06 

000007 

000008 

PHA 

PHA 

PHA 

PHA 

PHA 

PHA 

PHA 

PHA 

PHA 

PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 
PUT 


C1000  305 
L1010200 
CI000205 
C1J00205 
L 1030 200 
C1000205 
C100Q<t05 
L1020400 
C1000205 
L 1020 300 
C1000  305 
C1000405 
LI 040 300 
L10102Q0 
L 1020100 
LI  010 300 
L1030200 
L  1020  400 
L1020500 
L1040300 
L1Q50  400 

oej 

L1010200 
LI020 100 
L1J10200 
L1030200 
L1020400 
L1020500 
L 1040 200 
L1050  400 

C1000  201 
C1000201 
C10Q0  401 
C100050I 
C1J00102 
C1000  202 
C1000402 
C1000502 
C1000103 
C1000203 
C1000403 
C1000503 
C1J00104 
C1000204 
C1J00304 
C1J00504 
C1J00105 
C1000205 
C 1000 305 
C1000405 


■0.1000QE 
•0.100005 
•0.10000c 
•0.  10000E 
•0.10000E 
•0.10000E 
■0.10000E 
•0  .10000E 

0 
0 
0 

0 
0 

0 
0 
0 
0 
5 
0 
0 
2 
0 
0 
0 
0 
0 
0 
0 


02 
02 
02 
02 
02 
02 
02 
02 
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\PPENDIX  E 


OPTIMIZATION  PROGRAM,  INPUT  DATA 
MATRIX  PICTURE  AND  OUTPUT  FOR  MAX-SLACK  APPROACH 


//CHAN 

//* 

//* 

//* 

//   EX 

//ZPZ* 


2ND    RUN    FOR    MAX-SLACK1  ) 


G    JOB       (  1538,  1808 ),' CHANG    JOON    HOONG »  ,CLASS  =  A 

MPSIII     EXAMPLE    PROBLEM 

EC       MSSMPS 
SYSIN      00      * 

PROGRAM    <•  NO'  ) 

INITIALZ 

TITLE     {'EXAMPLE    PROBLEM, 

MOVE       (XOBJ,'QBJ'  ) 

MOVE        (XRHSi  •  INPUT1  J 

MOVE        UDATA, 'LI^EQS' ) 

MOVE        UPBNAME,  'PROJECT'  1 

CONVERT    <•  SUMMARY'  I 

SETJP     {'MAX') 

3CD0UT 

PICTURE 

WHIZARC 

PRIMAL 

SOLJTION 

EXIT 

PEND 


/* 

//EXEC.  SYSIN 
NAME 

ROwS 
N   OB 


CI 
CI 
CI 
CI 
CI 
CI 
CI 
CI 
CI 
CI 
Ci 
CI 
CI 
CI 
CI 
CI 
CI 
CI 
CI 
CI 
LI 
LI 
LI 
LI 
LI 
LI 
LI 
LI 
COLUMN 
XI 
XI 
XI 
XI 
XI 
XI 
XI 


E 
E 
E 
E 
E 
E 
z 
E 
E 

E 

c 

E 

E 

E 

E 

E 

E 

E 

E 

E 

E 

E 
c 

c 

E 
E 
E 
E 
E 


J 

0002  01 

000301 

000401 

000501 

000102 

000302 

000402 

000502 

000103 

0002  03 

000403 

0005  03 

000104 

000204 

000304 

000504 

000105 

0002  05 

000305 

000405 

010200 

020100 

010300 

030200 

020400 

020500 

0403  00 

050400 

S 

020101 

020101 

030201 

030201 

030201 

020401 

020401 


DO       * 
LINEQS 


C1J00201 
LI  020  100 
C 1000 201 
C1000  301 
L1030200 
C1000201 
C 1000 401 
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X1020401  L1020400 

X1020501  C1000201 

XI020501  C1J00501 

X1Q205Q1  L1020500 

X1040301  C1000301 

X1Q40301  C 1J00  401 

X1Q*0301  L1J40300 

X1050401  C1Q004Q1 

X1050401  C1000501 

X1050401  L1050400 

X101Q202  C10Q01Q2 

X1010202  L1010200 

X1010302  C1000102 

X1010302  C1000302 

X1010302  L1010300 

X1030202  C1000302 

XI 0302 02  L103O203 

X1040302  C10003O2 

X1040302  C10Q0402 

XI 0403 02  L1040  300 

XI 050402  C1000402 

X1050402  C10OQ5O2 

X1050402  L1050400 

X1010203  C10J0103 

X1010203  C1000203 

X1010203  L1010200 

X1020103  C1000  103 

X1020103  C1000203 

X1020103  L1020100 

X1010303  C1J00103 

X1010303  L1010300 

X1020403  C1000203 

X1020403  C1000403 

X1020403  L1J20400 

X1020503  C1000203 

X1020503  C1000503 

X1020503  L1020500 

X1040303  C1300403 

X1040303  L1040300 

X1050403  C10QQ403 

X1050403  C1000503 

X1050403  L10504Q0 

X1010204  C1000104 

X1010204  CI000204 

X1010204  L1010200 

X1020104  CI000104 

X1020104  C1000204 

X1020104  L1J20100 

X1010304  C1000104 

X1010304  C1000304 

X101Q304  L1J10300 

X1030204  C1000204 

X1030204  C1300304 

X1Q30204  L1030200 

X1020404  C1000204 

X1020404  Li02Q4Q0 

X1020504  C1000204 

X1020504  C1000504 

X1020504  L10205OO 

X1050404  C1J00504 

X1050404  11350400 

X101Q235  C1000105 

X10J.0205  C10G0205 

X1010205  L1010200 

X1020105  C1J00105 

X1020105  C1300205 

X1020105  L1020  100 


66 


X1010305 

01000105 

X1010305 

C1000  205 

—  1 

X1010305 

L1010300 

XI 0302 0  5 

01000205 

—  1 

XI 0302 0  5 

Ci 000 305 

X1030205 

L1030200 

X1020405 

01000205 

X1020405 

C1000405 

—1 

X1020405 

L1020400 

X1020505 

01000  205 

XI 0205 05 

L1020500 

X1040305 

01000  205 

—X 

X1040305 

01 000 405 

XI 0403 05 

L1040200 

S1000001 

OBJ 

L1010200 

1 

S1000002 

OBJ 

L1020100 

1 

SI  0000 03 

GBJ 

L1010300 

1 

S1000004 

OBJ 

L1030200 

1 

S1000005 

OBJ 

L1020400 

1 

S1000006 

OBJ 

L1020500 

1 

S1000007 

OBJ 

L1040300 

1 

S1000303 

OBJ 

L1050400 

i 

ALPHA 

OBJ 

* 

ALPHA 

L1010200 

-0, 

,  10000E 

02 

* 

ALPHA 

L1020  100 

-0, 

,  10000E 

02 

* 

ALPHA 

L1010300 

-0 

.10000E 

02 

* 

ALPHA 

LI  03  0  200 

-0. 

►10000E 

02 

* 

ALPHA 

L1020400 

-0 

.10000E 

02 

* 

ALPHA 

L1020500 

-0 

.10000E 

02 

* 

ALPHA 

L1040  300 

-Q< 

, 10000E 

02 

* 

ALPHA 

L1050  400 

-0, 

.10000E 

02 

RHS 

INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 

01000  201 
01000  301 
01000401 
01000501 
C1000  102 
01000202 
C1000402 
01000  502 
01000  103 
01000  203 
01000403 
01000  503 
01000  104 
01000204 
01000204 
01 000 504 
01000  105 
01000205 
01000205 
01000  405 
L1010200 
L1020  100 
L1010300 
L1030200 
L 1020 400 
L1020  500 
L1040  300 
L 1050 400 

0 
0 
0 

0 
0 

0 

0 

0 

0 

5 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

ENDA" 

/* 

// 
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APPENDIX   F 

LINK   FLOW    FOR  9/36  NETWORK 
Table   F.l.      Link   Flow   for  9/36  Network      (       );    MAX-SLACK 

,  TM„  DESTINATION  NODE 

LINK TOTAL 

(NODES)        12  3  4  5  5  7  8  9  FLOW 


1-2 

0.5 

(4) 
3.5 

(4) 
4 

2-1 

1-3 

(4) 

-t 

(4) 
4 

3-1 

1-8 

(4) 
3.75 

(4) 
3.75 

8-1 

0.5 

0.5 

1-9 

(3) 
3.75 

(3) 
3.75 

9-1 

2-3 

1.63 

1  .63 

3-2 

(1) 
2.25 

(1) 
2.25 

2-4 

(4) 
1  .88 

(4) 
1.88 

4-2 

(4) 
2.25 

(4) 
2.25 

2-9 

9-2 

(4) 

4 

(4) 
4 

3-4 

1  .88 

1.88 

4-3 

1.75 

1  .75 

3-5 

(4) 
3.75 

(4) 
3.75 

5-3 

(1) 
0.5 

(1) 
0.5 

4-5 

(4) 
3.75 

(4) 
3.75 

5-4 

4-6 
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i  TMI,  DESTINATION  NODE 

LINK       TOTAL 

(NODES)        12  3  4  5  6  7  8  9  FLOW 


6-4 

(4) 
4 

(4) 
4 

5-6 

6-5 

(4) 
4 

(4) 
4 

5-7 

7-5 

(D 
0.5 

(3) 
3.5 

(4) 
4 

6-7 

7-6 

(4) 
4 

(4) 
4 

6-8 

8-6 

(4) 
4 

(4) 

4 

7-8 

2.25 

2.25 

8-7 

1  .63 

1.63 

7-9 

(4) 
2.25 

(4) 
2.25 

9-7 

(3) 
1  .88 

(3) 
1  .88 

8-9 

1  .75 

1.75 

9-8 

1.88 

1.88 

REQ-MAT 

^(5)  = 

=15,    R5(2)=9 
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